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Abstract 

We compute the CP violation in the decays of heavy electroweak singlet neu- 
trinos, arising from both the one-loop vertex corrections and the wave function 
mixing. We extend the computation to the supersymmetric version of the model 
and discuss the implications for the generation of a lepton number asymmetry 
by the out of equilibrium decay of the heavy (s)neutrinos in the early Universe, 
to be reprocessed later in the observed baryon excess by anomalous electroweak 
processes. 



A very attractive scenario, proposed by Fukugita and Yanagida 0], for the genera- 
tion of the baryon asymmetry of the Universe, is based on the production of a lepton 
asymmetry by the out of equihbrium decays of heavy (m ^ TeV) electroweak singlet 
neutrinos. Their decays into hght leptons and Higgs bosons can violate CP if the Yukawa 
couplings involved have unremovable phases, and can then lead to the production of an 
excess of antileptons over leptons in the final state. The lepton asymmetry so produced 
is then partially converted into a baryon asymmetry by anomalous electroweak processes 
1^, which are in equilibrium at temperatures larger than the electroweak phase transition 
one, and as a result, the amount of baryons present in the Universe, ra^/ra^ ^ 5 x 10^^°, 
can be accounted for [|, ^. 

To postulate the existence of right-handed 'sterile' neutrinos constitutes one of the 
simplest and most economical extensions of the standard model, being strongly motivated 
by Grand Unified theories (GUT) such as SO(IO). It also allows to implement the see- 
saw mechanism which naturally accounts for non-zero, but still very small, light neutrino 
masses. Moreover, the values of the light neutrino masses suggested by the MSW solution 
to the solar neutrino problem, the atmospheric neutrino anomaly and the hot dark matter 
scenarios 0, point towards an intermediate scale, M ~ 10^-10^'^ GeV, for the right- 
handed neutrino masses. This further encourages the consideration of scenarios where 
the heavy neutrino decays act as sources for the generation of the baryon asymmetry 
of the Universe, since for these large masses it becomes easier to achieve the out of 
equilibrium conditions required for the efficient generation of an asymmetry, and also 
because temperatures comparable to the intermediate scale are more likely to be produced 
as a result of the reheating process at the end of inflation than the temperatures required 
in the conventional baryogenesis GUT scenarios. 

In addition to the non-supersymmetric version originally considered by Fukugita and 
Yanagida, the extension of this scenario to the supersymmetric case has been studied by 
Campbell, Davidson and Olive p, and also a related scenario where the asymmetry is 
produced in the decay of heavy scalar neutrinos produced non-thermally by the coherent 
oscillations of the scalar field at the end of infiation has been discussed by Murayama et 
al. 10]. A key ingredient for all these scenarios is the one-loop CP violating asymmetry 
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involved in the heavy (s) neutrino decay, and the reconsideration of this quantity will 
be the main issue of this work. As will be discussed, there are contributions to the 
asymmetry that have not been included previously, and we also want to discuss the 
different results present in the literature on this subject. 

Let us start with the non-supersymmetric version of the model, with a Lagrangian 
given by 



II = -Kjt^pN.PdtH^ + h.c. , (1) 

where = {ui l^) and = {H^ H^) are the lepton and Higgs doublets {eap = —^fSa, 
with ei2 = +1), and we are taking the Nj to be the Majorana mass eigenstate fields, of 
mass Mj. Since the scale of the heavy neutrinos is so much larger than the electroweak 
scale, it is reasonable to work directly in the symmetric phase where all particles except 
the heavy Majorana neutrinos Nj are massless, and all components of the neutral and 
charged complex Higgs fields are physical. 

The basic quantity we are interested in is 

I- Ni + I- Ne 

where Tm = Ea,/3 r(A^ Tif^) and T^^- = Ea,/3 r(A^ ^ i'^H^^) are the N decay rates 
(in the rest frame), summed over the neutral and charged leptons (and Higgs fields) 
which appear as final states in the N decaysQ. Note that at tree level 

Tiv.. = r^,, = ^^M,. (3) 

The asymmetry arises from the interference of the one-loop diagrams depicted 
in Fig. 1 with the tree level coupling. The vertex correction in Fig. la is the one that 
is usually considered, but it has however been pointed out p, ^, |10[ that the wave 
function piece in Fig. lb also contributes to the asymmetry, in an amount which is 
typically comparable to the vertex contribution. We have computed these CP violating 
asymmetries and we obtain for the vertex contribution to 

^For the Majorana light neutrinos, one should think of £ and £ as being the left and right-handed 
helicities of z/. 
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Figure 1: Diagrams contributing to the vertex (Fig. la) and wave function (Fig. lb) CP 
violation in the heavy singlet neutrino decay. 
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where Uk = M^/Mf, f{x) = y/x (1 — (1 + x)ln [(1 + x)/x]), and we have defined 



Im 



(A.A) 



(AtA), 



For the wave function piece we find 

Im{[M,(AtA)^^ + M,(AtA)jA*,A,4 



e^' (wave) 



(AtA), 



A useful check of the result can be obtained using the unitarity relation 



IT,; 



-2 Im 



if 



(5) 



(6) 



where we are retaining only the leading order terms; i and / are the initial and final 
states, and n are the possible intermediate states connecting them. T/„ are the transition 
amplitudes between the intermediate and final states, which in the present case arise 
from both the s and t-channel exchanges, which correspond to wave and vertex 
contributions respectively. 

The result in Eq. (4) coincides with the one recently obtained in ref. [0], differing by 
factors of 2 and 4 from those in ref. [11] and respectively^. The wave function piece. 



^The normalization of the e parameters are sometimes different (and not always explicit). 



3 



which in this scenario was only considered previously by Liu and Segre [0, is a factor 
of two larger than their result due to the fact that their computation applies to the case 
in which the scalar field is real and only the neutral lepton runs in the loop. Adding the 
charged lepton loop contribution the result is doubled. In Eq. (6) we do not include in 
the sum over the flavour of the intermediate neutrinos the case k = i, since in general 
states degenerate with the initial one do not contribute to the CP violating asymmetry 



T^. We also assumed that |Mfc — Mj| ^ T^^ in the computation. Notice that the 



contribution from k = i in the vertex piece, being proportional to Im {X^X)f- = 0, also 
vanishes. Hence, both sums in Eqs. (4) and (6) may be restricted to k ^ i. In Eq. (6) 
we have not summed over the final lepton flavour, but after summing over it one gets 

Q (^-a^^e) = --g^^^^-^J.., (7) 

so that the piece proportional to Mj in the square brackets in Eq. (6) actually gives a 
null contribution to the total lepton number asymmetry (although it may generate an 
asymmetry in the individual leptonic numbers Lj). In ref. ||10[ it was pointed out that 
there should also be in the vertex contribution, in addition to the result in Eq. (4), a 
term like the one just discussed appearing in Eq. (6). This term appears however only 
when a single real scalar Higgs field is allowed to run in the loop, but it can be shown 
that the new contributions arising from the two real scalars in the decomposition of the 
standard model Higgs fields if" = {H^ + iiff)/-\/2 actually cancel each other, leaving 
only the term proportional to Xki present in Eq. (4). 

In the particular case in which a hierarchy among the right-handed neutrino masses 
is considered, i.e. for ^ 1, we have that the wave function contribution becomes twice 
as large as the vertex one, and hence the total asymmetry produced in the decay of the 
lightest heavy neutrino^, Ni, becomes 

ef ^ {vertex) + {wave) ^ - A ^ ^ j,^. (g) 

The inclusion of the wave function piece then increases by a factor of three the CP 
violating asymmetry in the case of large mass hierarchies [|T0|. 

■^In the case of large hierarchies, the asymmetries produced by the heavier neutrino decays are usually 
erased before the lightest one decays. 



4 



If the departure from equilibrium in A^i decay is large, the lepton number asymmetry 
produced, per unit entropy, is 



S S Tl^ 



(9) 



where qn = "2 are the spin degrees of freedom of the Majorana neutrino A^i, so that their 
number density (assuming Maxwell-Boltzmann statistics) before they decay is gj^T^ /t^'^ 
if we assume that they were in chemical equilibrium before becoming non-relativistic (see 
below). Using that s = (2/45)5'* tt^T^, where (7* = gtosons + iy /^)gfermions is the effective 
number of relativistic degrees of freedom contributing to the entropy, {g.^ = 106.75 in the 
standard non-supersymmetric version of the model), we get 

^~4xl0-%^\ (10) 
s 

This lepton asymmetry will then be reprocessed by anomalous electroweak processes. 



leading to a baryon asymmetry ||13 



/24 + 4iVH\ 28 

Here is the number of Higgs doublets {Nh = 1 in the standard scenario considered 
now, while = 2 in the supersymmetric version to be discussed below if the scalar 
Higgs bosons are assumed to be lighter than the electroweak symmetry breaking scale). 
Combining Eqs. (10) and (11) and assuming that the Universe expansion is adiabatic, 
the present baryon asymmetry which results is 

^~-1.5xl0-='ef . (12) 

s 

Values of ~ —5 x 10~® are then required to account for the value ub/s c::^ 0.6-1 x 10^^° 
inferred from the successful theory of primordial nucleosynthesis | n | . In the case in which 
there is a hierarchy in the heavy neutrino masses, the Yukawa parameters need then to 
satisfy 

E ^ ^ki ^ 0.9 X 10-^ (13) 



Let us also briefly mention that the lepton asymmetry produced by decays is smaller 
than the value in Eq. (9) if the departure from equilibrium is not large during the decay 
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epoch, as is the case if the decay rate Fi = Tn^^ + Tj^j is comparable or larger than the 
expansion rate of the Universe when A^i becomes non relativistic, i.e. Fi > H{T = Mi), 
with H{T) being Hubble's constant at temperature T. We also want to emphasize that a 
problem of these scenarios is that the Yukawa interactions are not effective in establishing 
an equilibrium population of heavy neutrinos. To see this, note that before A^^i becomes 
non-relativistic, i.e. for x = Mi/T <^ 1, the thermally averaged rate scales as T~^, since 
(Fi) = Ki{x)/K2{x)ri ~ (a;/2)Fi oc Mf/T, where 7^1,2 are Bessel functions If. Taking 
into account that H oc T'^/Mpianck and supposing that Fi ^ H{T = Mi) (so that the 
departure from equilibrium during the decay is significant), one concludes that (Fi) ^ H 
holds true also for all temperatures larger than Mi. Hence, the inverse decays are unable 
to establish chemical equilibrium for A'^i at any temperature T > Mi. The pair production 
of Ni from Higgs boson or light lepton scattering may do somewhat better than inverse 
decays at high temperatures, since the relevant rate scales as {av) oc T, but it is anyhow 
insufficient to achieve chemical equilibrium for A'^i because the rates need to be quite 
suppressed at T ~ Mi in order not to erase any lepton asymmetry generated during the 
decay, and consequently are also suppressed with respect to the expansion rate at higher 
temperatures. Hence, additional interactions of the heavy neutrinos are required for 
them to have a thermal distribution prior to the decay, so that the leptogenesis scenario 
can be successful. This point was recently remarked in ref. [^], where it was shown that 
gauge interactions which are naturally present in GUT scenarios (new Z' or SU(2)/j gauge 
bosons) can easily do the job of producing a thermal population of heavy neutrinos at 
T > Ml. 

Turning now to the supersymmetric version of this scenario, the interactions of the 
heavy (s) neutrino field can be derived from the superpotential 

W = ^M,N,N, + \^e^pL^H^N,. (14) 

Supersymmetry breaking terms are of no relevance for the mechanism of lepton number 
generation and we are then left with the following trilinear couplings in the Lagrangian, 
in terms of four component spinors, 

C = -Xijeap S^MjNJLtH" + W^Pd^H^ + WYPJ^N, + (Pj^PiiV,L^| + h.c. . (15) 
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Figure 2: Supersymmetric contribution to the diagrams leading to CP violation in the 
N ^iH decay. 

From these couplings we obtain the tree level relations 

r^v^^ + r^^^- = r^^^ + r^^^. = r^.^ = r^^ = ^—^^Mi. (i6) 

There are now new diagrams, like the ones in Figure 2, contributing to the generation 
of a leptonic asymmetry. We will denote the corresponding asymmetry parameters in the 
supersymmetric case with a tilde, so that for instance (vertex) will receive contribu- 
tions from the 'standard' diagram in Figure la as well as from the one in Figure 2a which 
involves superpartners in the loopQ. We also define the slepton asymmetry associated to 
N decays as 

i - r - _L r - ' ^ ' 

^ NL^ ^ NL* 



which arises from diagrams similar to those in Figures 1 and 2 but with Lh in the final 
state. We similarly define the asymmetry parameters associated to sneutrino decays as 



~N* _ ^ N*l ^ Nl ~N _ ^NL ^ff*L* /i t 

^ ~ F- +r-- ' ~ r- - + r- - ■ ^ 



After a direct computation of these quantities we obtain 

1 



' (vertex) = -—J^d iVk) ^ki , (19) 
°^ k 



^For later convenience, we define as in Eq. (2), rather than normahzing it to the total N decay 
rate which also includes the slepton final states. 
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where g{x) = a/x In [(1 + x)/x], while 

i^{wave) = 2ef{wave). (20) 



For the remaining channels we get similar results 

(vertex) = (vertex) = (vertex) = (vertex), (21) 
and also the wave function contributions are equal 

(wave) = (wave) = (wave) = (wave). (22) 

The vertex pieces were computed previously in ref. and our results are propor- 
tional. The sneutrino decay asymmetry was also considered in ref. in the limit in 
which all heavy (s) neutrinos are degenerate and, specialised to that case, our results are 
in agreement (except for the overall sign). The wave function contributions were not 
considered previously and, as in the non-supersymmetric case, they are non-negligible. 
Let us also note that there are additional one loop diagrams involving the four-scalar 
couplings appearing in the F-terms of the scalar potential. Although helpful to cure the 
divergences in the loops, they do not contribute to the asymmetry generation. 

For the study of the decays of thermal populations of heavy neutrinos and sneutrinos, 
it is convenient to introduce the total asymmetry defined as 



+ef +ef 



(23) 



resulting in 



--E 



givk) + 



I, 



ki ■ 



(24) 



In particular, in the case of hierarchical masses (yu ^ 1), we find that again the wave 
contribution becomes twice as large as the vertex one, giving 



3^1 

kytl 



27T — ^/yk 



(25) 
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With the definition in Eq. (23), the resulting lepton asymmetry i^ 

^ = ^ 1 X 10-%, (26) 

S S Tl'^ 

where we have used that the effective number of degrees of freedom in the supersymmetric 
scenario is approximately doubled, i.e. gf^^^ ~ 2 gf^'. Hence, to account for the baryon 
asymmetry of the Universe we need now 

E ^ ^ 0.7 X 10-^ (27) 

Comparing this with Eq. (13), we see that the required parameters are similar in the 
supersymmetric and standard scenarios. 

To summarize, we have computed all the contributions to the CP violating asym- 
metries arising at one-loop in the decays of heavy (s)neutrinos, both in the standard 
non-supersymmetric Fukugita-Yanagida scenario and in its supersymmetric version. We 
have discussed the different results present in the literature and showed that the con- 
tribution from wave function mixing is relevant in the computation of the CP violating 
asymmetries. The baryon number generated in both scenarios was also obtained. 

It is a pleasure to thank A. Masiero for helpful suggestions and S. Bertolini, N. Paver 
and A. Smirnov for discussions. 



^Had we normalised to the total N decay rate, their contribution to would have to be 
multiplied by (7Ar = 2 as in Eq. (9). 
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